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Abstract. We study the relation between entanglement and quantum phase 
transition (QPT) from a new perspective. Motivated by one's intuition: QPT is 
characterized by the change of the ground-state structure, while entangled states 
belong to different classes have different structures, we conjecture that QPT occurs 
as the class of ground-state entanglement changes and prove it in XX model. Despite 
the classification of multipartite entanglement is yet unresolved, we proposed a new 
method to judge whether two many-body states belong to the same entanglement class. 
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1. Introduction 

Quantum phase transition (QPT) is a phase transition that occurs at absolute zero 
temperature, and means nonanalyticity of the ground-state properties. The singularity 
may be a discontinuity in the first or higher-order derivative of the ground-state energy, 
respectively referred as first-order or continuous QPTpQ. 

QPT is usually accompanied by a qualitative change in the nature of the correlation 
in the ground state, and quantum fluctuation is the ultimate reason leading to QPT, 
so it is certainly of major interest in both condensed matter physics and quantum 
information science to describe the connection between QPT and quantum entanglement 
[2]. Various entanglement measures are calculated and hoped to exhibit singular 
behavior at quantum critical point. Some bipartite entanglement measurements such 
as concurrence (3j IU [5], entanglement entropy j6j [7], can indeed identify particular QPT. 
There even a general theory of the relation between QPT and bipartite entanglement 
was developed under certain conditions [E]. However, counterexample was soon foundpj. 
Therefore, QPT in terms of multipartite entanglement began drawing attentions |10[ ITT] . 
It is also a noteworthy problem whether multipartite or bipartite entanglement being 
favored at QPTpJ Q3]. 

So far there is not a universal conclusion of the relation between entanglement 
(bipartite or multipartite) and QPT (first-order or continuous). Considering that the 
intrinsic character of QPT is the change of the ground-state structure, so the problem is 
probably due to that the singularities of entanglement measures aren't essential to the 
change of the ground-state structure. Geometric phase [Hj or fidelity [151 EH EH] ma Y 
be a good indication of QPT. But concerning entanglement, we think that inequivalent 
entanglement class can appropriately reflect the change of ground-state structure, since 
two inequivalently entangled states have different structures. So a promising way to 
reveal the deep connection between QPT and entanglement is to study the classification 
of ground-state entanglement around QPT. Our focus in the paper is to investigate the 
classification of entangled ground states in the vicinity of QPT in one- dimensional XX 
model [18]. It is shown that, no matter the model length N is arbitrary or tends to 
infinity, the change of ground-state entanglement class always indicates the occurrence 
of QPT. 

2. SLOCC classification of entanglement 

As concerning the classification of entangled states, stochastic local operation and 
classical communications (SLOCC) are usually used to define equivalent classes. That 
is, two states are said to belong to the same entanglement class if both of them can be 
obtained from the other with nonzero probability by means of local operation assisted by 
classical communications. Many researchers have investigated the SLOCC-inequivalent 
classes of pure entangled states [191 12Q1 ISJJ . 

The complete classification of pure entangled states is indeed an intricate task. 



Classification of entanglement and quantum phase transition in XX model 



3 



However, if we only need to judge whether two states are SLOCC-inequivalent, there is 
a simple criteria by virtue of Schmidt decomposition [2]. As well known, for any bipartite 
pure state G H = Ha ®1~Lb, there exist local orthonormal bases {!««)} G %a and 
{\vi)} £%b such that 

\9) = Y i a i \u i )®\vi), (1) 

i 

a; are nonnegative real numbers satisfying YU a f — 1, referred to as Schmidt coefficients. 
The number of nonzero a, is known as the Schmidt rank, here denoted by Sch(|\l/)). 
It can be easily deduced that any two SLOCC-equivalent states must have the same 
Schmidt rank j22j E3]. In other words, two bipartite states with different Schmidt 
rank are SLOCC-inequivalent. It provides a clue to judge the SLOCC-inequivalence of 
multipartite state. Given two iV-party states, we can calculate its Schmidt rank based 
on a particular bipartition. A bipartition means a division of iV-party system into two 
nonempty and disjoint parts, i.e., one part including M(l < M < N) bodies and the 
other (N — M) bodies. If two iV-party states have different Schmidt rank based on the 
same bipartition, they necessarily belong to different entanglement class. This method 
bypasses the involved issue of complete classification of multipartite states and may be 
crucial in the researches of QPT in many-body system. 



3. First-order and continuous QPT in XX model 

The Hamiltonian of XX model is 

j N N 

u = j EK^ 1 + 4< +1 ) ~ B T.<- (2) 
4 i=i i=i 

where o%/ y / z are the usual Pauli matrices of the zth spin (cyclic boundary condition 
N + 1 = 1 is assumed). The external magnetic field B could always be supposed 
positive without loss of generality. The model can be analytically solved by Jordan- 
Wigner transformation [18]. Using the operators 
i fc-1 
a± = -(a x ± io-y), 4 = <r* IJ er*, 

/ 8=1 

the Hamiltonian (|2J) is transformed into 

J JV ~ 1 J N 

U = -- + c Ui+i) + -a(c\c N + 4d) -BN + 2BJ2 4^ 

1 i=i z i=i 

where a = n&Li(l — ^ c k c k) — Uk=i a z = ( — l) r - r * s the total number of spin-downs 
which is a constant. Introducing the Fourier transformation of c\: 

v iv k=1 Jy 
where n is integer (half-odd integer) for odd (even) r, the Hamiltonian is that of one- 
dimensional spinless ferminons 

U = -BN + J2(2B - J cos q)C^C q . 
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The lowest energy eigenstate with fixed r can be expressed as 



r 



lV>oO=n^(r + l-20/iv!£l^}> 



(4) 



1=1 N 



and its energy is 



r 



El = - Jj^ cos[7r(r + 1 - 21) /N] - B(N - 2r) 



i=i 




= -D r J-B(N-2r). 

Obviously,-D r = D r . Thus when J is fixed and B is tuned, 

i N 



(5) 



includes all of the possible ground states of the XX model. With the modulation of 
the external magnetic field, each state in Q becomes the ground state in turn, while the 
energy of the system changes abruptly. The first-order derivative of the energy with the 
magnetic field B is discontinuous. A first-order QPT occurs. 

For iV-limited XX model, the first-order QPT occur at critical values of 
magnetic field, which can be achieved from the eigenenergy corresponding to each state 
in Q. That is, the energy of \i\f^) is —D r J — B(N — 2r), then the transition from |^q) 
to l^o 44 ) occurs at 



For iV — > oo, it is already known that at B c = J/2, there occurs a continuous QPT, 
i.e., a superfluid-Mott insulator phase transition pQ. 

4. SLOCC classification of ground-state around QPT in XX model 

Now we begin to investigate the ground-state entanglement class around QPT in XX 
model. 

When N — > oo, there is only one contimuous QPT. At B c = J/2, the system 
transits from Mott-insulator phase to superfluid phase or vice versa. The ground state 
of Mott-insulator phase is a separable state with all spins pointing to the same direction, 
while the ground state of superfluid phase is sure to be an entangled state [Tj. Therefore, 
the class of ground-state entanglement changes when QPT happens. 

When N is limited, since every stats in Q could be ground state with the adjustment 
of magnetic field, we need to prove that all states in Q is inequivalently entangled. 

To complete this tough task for arbitrary N, we calculate the Schmidt rank of every 
element in Q, based on [yj <8> N — [yj bipartition. We will show that Sch(|V>o)) — 2 r 
holds for arbitrary r and N. For simplicity, [yj =M henceforth. 



-D r J + 2Br 



D r+1 J + 2B(r + l), 



i.e., 




r = 
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Notice that 



44|o_-_o> 

N 



\mn), 
-\mn), 
0, 



if m < n, 
if m > n, 
if m = n, 



where \mn) corresponds to spin configuration in which all spins are up, except the spin 
at the site m and n are down. Replacing Eq.(j3]) into (jlj), the ground state with r 
spin-downs can be expressed as 



r N 



IIE ex p 

1=1 k=l 



i(r + l -2l)kn 



N 



4 1 0^0) 

N 



■ — exp 

Vjv i<fe 1 <...<fe r .<jv 

. > sin 

- / A/7 LL 

V^V Kt 1 <...<fc r <JVKi<j<r 



i(r — l)kin 




i(r — 3)/c27r 




i(l — r)k r n 


exp 


■ ■ ■ exp 


N 


N 


N 



(ki - kj)TT 



N 



| k\ ' ' ' . 



I** 



The last step is achieved using e lx — e 1X = 2i sin x. 



Next we calculate the Schmidt rank of based on M ® iV - M bipartition. The 
constant (2i) c ^/y/N^ can be omitted. \i)jq) is a weighted superposition of all possible 
\k\ - ■ -k r ) where the value range of all ki are [1, N] and ki < • • ■ < k r must be satisfied 
, so we rewrite |^q) as 



M 



e n sin 

M<ki<...<k r <N l<i<j<r 



(ki - kj)lT 



N 



| k\ ' ' ' kf ) 



+ •■■+ E 

l<fci<---<fc ; <M 



\h---ki)<S> II sin 

M<k l+1 <---<k r <N l<i<i<r 



+ •■■ + 



e n sin 

l<fci<---<fc r <iW l<i<j<r 



(ki - kj)ir 



N 



(ki - fcj)7T 

iV 

0_^_0) 

N—M 



That is, we first carry out a preliminary Schmidt decomposition by divide all 
possible \ki - ■ ■ k r ) into r + 1 groups, according to the number of spin-downs that locate 
in the former M qubits, i.e., 



hAo> = EE4K>®K4 

1=0 i 



(6) 



Obviously, when k ^ I, Vi, j, (u^u 1 ) = (vf\v l j) = always hold. If Sch(\ipQ^)) represents 



the number of nonzero a;-, then 



Sch(|^)=£Sch(|^>). 

1=0 



We will explain that 

Sch(|Vo W » = C, 



(7) 



(8) 
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thereby 



Sch(|^» = £^ = 2 r . 



(9) 



1=0 



Notice that the maximum value of r is M, so the Schmidt rank of \ip r ) will never 
exceed 2 M 

First it is easily found Sch(|^o (0) » = Sch(|$J (r) )) = 1. Eq. © holds for l = 0,r. 

Next remembering that if {ai, - ■ ■ ,a{\ is a set of linear independent vectors, 
then we can achieve an equivalent set of orthogonal vectors {Pi, ■ ■ ■ , Pi} by Schmidt 
orthogonalization. Thereby, Sch(|^o^)) is the rank of such a C l M <8> (Tridimensional 
matrix A r<yl \ Every element of A r "> can be uniformly expressed as Tli<i<j<<r si n [(^i — 
kj)ir/N], i.e., a product of C 2 sine functions. For each row 1 < ki < ■ ■ ■ < ki < M are 
fixed and M < ki + i < ■ ■ ■ < k r < N vary, while for each column the situation is just 
the reverse. We find the rank of A r ^ by elementary row(column) transformation. So 
every element can first be simplified as rii<i<;,z+i<j<r sin[(fci — kj)ir/N], i.e., a product 
of l(r — I) sine functions, as for each row of A r( - l \ IIi<j<7<z sin[(/cj — kj)ir/N] is a 
constant, and for each column llz+i<i<j<r sin[(/cj — kj)n/N] is a constant. Furthermore, 
Sch(|^ (1) )) = Sch(|^o (r_1) )) should hold. The reason is that the elementary row and 
column transformation of A r ^ and A r ^ r ~ 1 ^ will yield similar simplest form. Then we 
only need to calculate Scb-d^o^}) by elementary row transformation for I — 1, • • • , |_§J ■ 

Take r = 2, 1 = 1 as an example, every element can be expressed as sin[(A;i— k-zjit/N], 
and for each column 1 < ki < M varies while for each row M < k% < N varies. So the 
the matrix is (overall minus sign is omitted) 



^2(1) 



sin 



N ' 



sin 



sin(f) 



'2tt\ 
, N I 

sin(f) 



sin(fTr) sin(M±i7r) 



sin(^Mvr) 



sin(%±7r] 



N 



Let the ith row vector of A 2 ^ is denoted because 



/ rrni 

sin Ur 

so for i = L...M-2, 



+ sin 



(m + 2)7r 
N 



7T 



2 cos ( — ) sin 



(m + 1)tt 
N 



(10) 



a* + ai+2 = 2 cos ( — | a,+ . 



7T 



always holds. Then by elementary row transformation a i+2 =aj +2 + a^ — 2 cos aj +1 
for i = 1, • • • , M — 2, can be transformed into 



A 2(D 



' Af-1 . 
> AT 

sm(fTr) 



sm 



AT 
M+l 

. AT 



7T 



sm(^vr) 
sm(^Tr) 
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So the rank of A 2 ^ is 2. Elementary column transformation apparently yields the same 
result. 

For bigger r and I, although the matrix A r ^ becomes complicated rapidly, the 
knack used to find the rank of A r ^ is analogous, however, much more intricate. We will 
expatiate step by step. 

First consider A r<yl \ every element is in the simplified form of FI 2 <j<r sin[(A: 1 — 
kj)7r/N], i.e., a product of r — 1 sine functions. It can be transformed into a 
linear sum of a series of sine functions like sin{[pA;i + f p (k2,- - • ,k r )]ir/N}, where 
p — r — 1, r — 3, • • • , 1(0) for even (odd) r. After elementary row transformation like 
above, there are two rows left for every nonzero p and only one row left when p = 0. 
Then the simplest expression of A r ^ is 



(M-r+l-kj)ir 
N 



sin 



(M-kj)iT 



N 



A r(D 

Il2<7<r 

so the rank of AW» = r = C}. 

Next consider A r<y2 \ the element is simplified as Y\ 3< j <r sin[(/ci — kj)%/N] sin[(A;2 — 
kj)ir/N], which is a product of 2(r — 2) sine functions. When we simplify A r ^ by similar 
strategy, we must bear in mind that both k\ and k 2 vary in the value range [1 , M] and 
k\ < &2 must be satisfied. So the simplest form of A r ^ is 

~ {M-r+l-kj)ir' 



A r(2) 



n 



3<J<r 



sin 



iY 



sin 



(M-r+2-kj)n 
N 



n 



3<i<r- 



sin 



(M-r+l-kj)iT 



Il3<j<r Srn 



N 



(M-l-kj)-w 



sm 



(M-kj)iT 



N 



N 



sm 



(M-kj)-K 



N 



the rank of A r ^ = 1 + 2 + • • ■ + r - 1 = C 2 T . 

Finally, analogous but involved generalization yields the ultimate expression of A r ^ 

is 



A r(D 



n 



l+l<j<r 



sm 



(M-r+l-kj)-K 



N 



■ ■ sm 



(M-r+l-kj)ir 



Il/+l<7<r S i n 



(M-l + l-kj)TT 



N 



■ sm 



N 



(M-kj)-K 



N 



and its rank is 



(r - / + 1) + (r - /)(1 + 2) + (r - I - 1)(1 + 2 + 3) 
+ --- + (l + 2 + ---r-/ + l) 



r E 1 [r-(/-2+ J )] 

Cl 



i=l 



So far we have explained Eq. ([8]), accordingly Eq.Q is proved for arbitrary r and 
N. So we have proved that in arbitrary iV-qubit XX model, no matter N is infinite or 
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finite, the occurrence of all QPT, first-order and continuous, can be witnessed by the 
change of ground-state entanglement class. 

5. Conclusion 

In conclusion, we have studied the relation between classification of ground-state 
entanglement and QPT in A-qubit XX model. For arbitrary A, when the exchange 
constant J is kept invariable and the external magnetic field B is tuned, QPT occurs. 
For N — > oo or limited N, the QPT is continuous or first-order respectively. No matter 
what the type of QPT is, we find that the occurrence of QPT is always indicated by 
the change of class of entangled ground state. Although the conclusion obtained in the 
XX model seems too particular, we think it is indeed a reasonable conclusion that 
the entangled ground states in the vicinity of the transition are SLOCC-inequivalent. 
Because the intrinsic feature of QPT is the change of the structure of the ground state, 
and inequivalently entangled states have different structure. We believe our results grasp 
the essence of the relevance of entanglement in QPT and hope it can be verified broadly 
in the future. Besides, the proof based on the Schmidt rank provides a partial solution 
to judge SLOCC-inequivalent entanglement. Although it works only if Schmidt rank 
indeed changes at the QPT, since the coincidence of the Schmidt rank does not ensure 
the same entanglement class for three or more components, it develops a new method 
regardlessly the complete classification of multipartite entangled is far to be resolved 
nowadays. 
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